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Abstract 

The goal of this work is to determine classes of traveling solitary wave 
solutions for a differential approximation of a finite difference scheme by 
means of a hyperbolic ansatz. 



1 Introduction 

The Burgers equation: 



u t + cuu x - fiU xx = 0, 



a, fi being real constants, plays a crucial role in the history of wave equations. It 
was named after its use by Burgers [I] for studying turbulence in 1939. 

A finite difference scheme for the Burgers equation can be written under the fol- 
lowing form: 

F(u l m ,h,r) = 0, (2) 

where: 

it™ = u (I dx , m dt) (3) 
* Corresponding author: david@lmm.jussieu.fr; fax number: (+33) 1.44.27.52.59. 
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I G {i — 1, i, i + m G {n — l,n,n+ 1}, j = 1, n x , n = 1, n 4 , h, r denot- 
ing respectively the mesh size and time step, and a the Courant-Friedrichs-Lewy 
number (cfl) coefficient, defined as a = cr/h. 



A numerical scheme is specified by selecting appropriate expression of the function 
F in equation (J2J. 

Considering the ui m terms as functions of the mesh size h and time step r, ex- 
panding them at a given order by means of their Taylor series, and neglecting the 
o(t p ) and o(h q ) terms, for given values of the integers p, q, lead to a differential 
approximation of the Burgers equation, of the form: 

d r u d s u . 

r, s being integers. 

For sake of simplicity, a non-dimensional form of (J3J will be used: 

^,p,p) = 0, (5) 
y dx r dt s 

Depending on this differential approximation (jlj), solutions, as solitary waves, may 
arise. 



The paper is organized as follows. Two specific schemes are exhibited in section 
El The general method is exposed in Sectional In Section |U it is shown that out 
of the two studied schemes, only one leads to solitary waves. A related class of 
traveling wave solutions of equation (jlj) is thus presented, by using a hyperbolic 
ansatz. The stability of this class of solutions is discussed in section A numerical 
example is exposed in section H3 



2 Scheme study 

2.1 Finite second-order centered scheme in space, 
Euler-time scheme 

For the finite second-order centered scheme in space and Euler-time scheme, the 
function F of (j2J) takes the form: 



w m. x Ui n+1 - Ui n n U i+l n -U i _ l n U i+l n -2u,r + U l . l n 

F(ui m ,h,r) = + c Ui n — n = (6) 
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Consider function of the time step r, and expand it at the second order 

by means of its Taylor series: 



r 2 

Ui n+1 — u (i h, (n + 1) t) = u {i h,nr) + r u t (ih,nr) H w tt {i h,nr) + o(r 2 ) (7) 



It ensures: 



Uj n+1 — Uj n . , . T 



«t (i/i,nr) + - (i h,nr) + o(r) 



r »v > ■ 2 

In the same way, consider u^+i™ and Wi_i n as functions of the mesh size h, and 
expand them at the fourth order by means of their Taylor series: 

iti+i™ =u ((i + 1) h, n t) 

h 2 h 3 h A ( 9 ) 

=u (ih,riT) + h u x (ih,nr) + — (i h, n r) + — u xxx (ih,nr) + — (i h,nr) + o(h 4 ) 



=m ((i — 1)Ii,iit) 

h 2 h 3 h 4 ( 10 ) 

=m (ih,riT) - h u x (ih,nr) + — {ih,nr) - — u xxx {ih,nr) + — ti IIra (j h,nr) + o(/i 4 ) 

It ensures: 

Mj+l™ — 2«i n + Mj-i™ 2h 2 2 

— = it a; a; {ih,nr) + —^-u xxxx (ih,nr) + o{h ) (11) 

and: 

l+ 2h = u x (ih,nr) + — u xxx (i h,nr) + o(h 3 ) (12) 
Equation (£Q) can thus be written as: 

T 

ut (i h,nr) + — utt (i h,nr) + o(r) 

h 2 

+ cu(ih,nr) [u x (ih,nr) + — ti IJ;I (i h,nr) + o(/i 3 ) ] (13) 
2 /i 2 

- [«a;a: (« fr, n "7") H — Mn^a; (« h,nr) + o(h 2 )] = 

i. e. , at x = i h and t = nr: 



[«t + ^ utt + o(t) + cu [u x + ^- + o(/i 3 )] - fi \u xx + u xxxx + o(h 2 )] ] = (14) 
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The first differential approximation of the Burgers equation (JH) is thus obtained 
neglecting the o(r) and o{h 2 ) terms: 

\ u t + - u t t + cu \u x + — u xxx \ - \u xx + — u xxxx ~\ 1 =0 (15) 
L 2 L 3! JL 12 JJ (®,t) 

that we will keep as: 

r h 2 h 2 

u t + cuu x - au xx + -u tt + —uu xxx - fx — u xxxx = (16) 

2 o 12 

For sake of simplicity, this latter equation can be adimensionalized through in the 

following way: 

set: 

u = UqU 

t = T t (17) 

x = ho x 



where: 



T 

Denote by Reh the mesh Reynolds number, defined as: 



U = ^ (18) 



Unh , . 

Re h = — (19) 



The change of variables (|T7jl leads to: 



(20) 



Multiplying JEJ by ^ yields: 



~ _l U ° T ° ~ ~ _ r ° ~ r ~ ^ U o T o ~ ~ _ fe2r o ~ 

u t ' c 7 uu x A 4 ,2 M ^ o M tt ~r „ ,3 It W^s^ 4 M 

Relations (JTBJ) and (|T9"j) ensure: 



h . r _ /i 2 fr 3 

~\~ CU U x U xx + It^r + 2 It ttjjaja; „ Uxxxx — U (22) 

riotieh 2 To o n llnehiio 

For h = fin, due to a = ^r 21 , (1231) becomes: 



1 1 „ _ 1 

Ul + CUU X — tt^j + (J Ull ~l~ 77 u u xxx 77T~n u xxxx = (23) 

ite^ 12 Reh 
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2.2 The Lax-Wendroff scheme 

For the Lax-Wendroff scheme, the function F of (0) takes the form: 

F(«, ,M = + cu , | )( ^ f" = 

(24) 

h 2 



— — is expressed by means of (jBJ), and — 2 ^ +M ' 1 by means of pijl : 
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u i+ i n — 2 Ui n + Ui-i n . 2h 2i 
— = u xx {i h,nr) + u xxxx (i h,nr) + o{h ) (25) 

fTTj) also yields: 

l+ 2h = u x (ih,nr) + — u xxx (i h,nr) + o{h 3 ) (26) 
Equation can thus be written as: 

T 

ut (i h,nr) + — utt \i ft., nr) + o(r) 

ft 2 

+ au(ih,nr) [u x (ih,nr) + — tim (i ft, n t) + o(ft 3 ) ] (27) 

c 2 t 2 ft 2 

— (jj, H — ) [w^ (i ft, nr) + — - (i ft, nr) + o(ft 2 )] = 



i. e., at x — i h and t = nr: 



I" u t + -u t t + o(r) + cm [uz + — u xxx + o(ft 3 )l — (/i H — ) \u xx + —— u xxxx + o(ft 2 )l 1 = (28) 

The first differential approximation of the Burgers equation (JIJ) is thus obtained 
neglecting the o(r) and o(/i 2 ) terms: 

[«i + -it(+cu[a I + — u„J — (/-H — ) \u xx + — iiuJ 1 = (29) 

that we will keep as: 

C 2 T T /i 2 C 2 T /l 2 

w t + cuu x - (fi + jj^) u xx + -u tt + —uu xxx - (/i + — ) —u xxxx = (30) 
Equation (J3Uj) is adimensionalized as in section 12.11 leading to: 

,1 c 2 o\ r ^ 1 „ „ .1 c 2 o\ 1 _ n /ni , 

Ut + cuu x - ( — h -— J u xx + - u~ti + - uu xxx - ( — h -—-) —u xxxx = (31) 

Keh 2 2 b Keh 2 12 



3 Solitary waves 



Approximated solutions of the Burgers equation (0) by means of the difference 
scheme (J2J) strongly depend on the values of the time and space steps. For specific 
values of r and h, equation (jHJ) can, for instance, have traveling wave solutions 
which can be of great disturbance to the searched solution. 

We presently aim at determining the conditions, depending on the values of the 
parameters r and h, which give birth to traveling wave solutions of (j!6|) . 
Following Feng [2] and our previous work [3] , where traveling wave solutions of the 
cBKDV equation were exhibited as combinations of bell-profile waves and kink- 
profile waves, we aim at determining traveling wave solutions of (JHJ). 
Following |2j, we assume that equation (jHJ) has the traveling wave solution of the 
form 

u(x,t)=u(£), £ = x-vt (32) 
where v is the wave velocity. Substituting (}32j) into equation (jSJ) leads to: 

^(u,u {r \(-v) s u {s) ) = 0, (33) 

Performing an integration of (}3*3*j) with respect to £ and setting the integration 
constant to zero leads to an equation of the form: 

^(u,u^,(-v) s u is) ) =0, (34) 
which will be the starting point for the determination of solitary waves solutions. 

4 Traveling Solitary Waves 
4.1 Hyperbolic Ansatz 

The discussion in the preceding section provides us useful information when we 
construct traveling solitary wave solutions for equation (jHHJ). Based on this result, 
in this section, a class of traveling wave solutions of the equivalent equation (|TB|) 
are searched as a combination of bell-profile waves and kink-profile waves of the 
form 

n 

u(x,t) = ^2(U i tanh i [d(x-vi)] +1^86^ [d(x - vt +x ) ]) + V (35) 
i=i 

where the U-s, V-s, Cjs, (i — 1, ■ • • , n), Vq and v are constants to be determined. 
In the following, c is taken equal to 1. 
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4.2 Theoretical study 

Substitution of (J33)) into equation (|3~3*)) leads to an equation of the form 



M tanh 4 (C, sech J (Q f ) sinh fc (Q f ) = (36) 
the Ai being real constants. 

The difficulty for solving equation (}36|) lies in finding the values of the constants 
Ui, Vi, Ci, Vo and v by using the over-determined algebraic equations. Follow- 
ing |2j, after balancing the higher-order derivative term and the leading nonlinear 
term, we deduce n — 1. Then, following [3] we replace sech(Ci£) by eCl g + 2 e -<j lS , 

sinh(Ci^) by eCl ~ e Cli , tanh(Ci £) by ^4^^4 > an d multiply both sides by 

(e^ 1 + e - ^ 1 ) 5 e 5 ^ 1 , so that equation (|36p can be rewritten in the following form: 

10 

J2 p k(Ui, V u d, v, V )e kC ^ = 0, (37) 

fc=0 

where the Pk (k = 0, 10), are polynomials of Ui, V±, Ci, V and v. 
4.3 Scheme study 

4.3.1 The Finite second-order centered scheme in space, Euler-time 
scheme 

Equation (|3*3|) is presently given by: 



-v «'(£) + c 5(0 u'(0 - 5"(0 + « 2 £ u"(0 + ~ 5(0 u {3) (0 - i ^ (4) (0 = 

(38) 

Performing an integration of (|38|) with respect to £ and setting the integration 
constant to zero yields: 



(39) 

The related system (|37j) has consistent solutions, which are given in Tables d 
For sake of simplicity, we use e to denote 1 or —1. 
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Table 1: 



Set 4 



o 

C„ f „ 1 o 484jfe h 108 108 Fj ~ 6 108 ~ 

° eLb Z 729 fc llv^H-Reft fc 5yiTiJe h U b Vn b 5v/IT Re h 

Set3 ^cl-y _ifL o Gl ^ 

6C 2 (4C 2 -9) 5ifeh 5fieh 5 Reh 

fe fc (64C«-384C*+551C?-156) "Tj 5 ^" 0,1 2(8Cf-9Ci) i 8Cl 

602(402-9)" ^ fle fe (8C2-13) J?e h (8C 2 -13) 



In the following, we shall denote: 

484 Re h 
729 

5Re h (l7C* 2 -12) 

6C2(4C2-9) 2 (40) 
ifc h (64C 1 6 -384C 1 4 +551C 1 2 -156) Re h (C 2 -4) (8C?-13)(8C?-3) 

6C 2 (4C 2 -9) 2 6C 2 (4C2-9) 2 

4.3.2 The Lax-Wendroff scheme 

Equation is then given by: 

-vu'(0 + c«({) u'(0 - ( J- + 4^) fi"(0 + fi"(0 + i fi(0 « (3) (0 - + 4^) 4 « (4) («) = ( 41 ) 

v Re^ 2 ' 2 6 Reh 2 12 

Performing an integration with respect to £ and setting the integration constant to 
zero yields: 

-v m + \ u\i) + (,2 1 _ ( _L + ^)) u ' (0 + I [ fi(£) fi w (0 _ I (0 ] _ (_L + J_ fi (3) (0 = o 

{42) 

The related system (pTTj) does not admit consistent solutions. 

5 Stability study 

In the following, the stability of the solutions presented in section l4*.3.1l is discussed. 

The variations of the cfl coefficient a as a function of the parameters Ret, Cx, 
have crucial influence on the stability. 



Cl,2 
0-3 

04 
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(73, o"4 being even functions of C\, we shall restrain our study to C\ G ] 0, +00 [. 
Calculation yield: 



20(34 C*i 4 -36C*i 2 +27) 0[<T 3 -_R.e h ] 20(34 d 4 -36 C 1 2 +2l) Re\ 

; 7^ ^ U , o775 = ; ^ U 

3C 1 3 (4C 1 2 -9) ' 8C i 3C 1 3 (4C 1 2 -9) 

4(96 Cf -238 C 4 +468 C* 2 -35l) Re h d\u 4 -R eh \ 4(96 Cf-238 C 4 +468 C 2 -35l) Re\ 

3C3(4C 2 -9) 3 ' ^C 2 ~ 3C 3 (4C 2 -9) 3 

(43) 

Denote by C® 2 the value of C\ for which ^ h - and - ^ ^ ^ vanish. 

For a given value of the mesh Reynolds number Reh, we obtain the following 
interesting variation tables: 



c\ 


f +oo 


d [R^\ 9 vs-Reh 


+ II + 


dC'f ' dC'f 


Rlk > ^ ■ 


+oo + oo 

/ \ 

-oo || 



(44) 



c\ 


Ci°" | +oo 




+ o - || + 


dC'f ' <9C* 2 




+oo +oo + oo 

\ / \ 

< < 



(45) 



a Reh take all the values between and +00. Thus, for stable and unstable 
sets (cr • Reh), there exists a solitary wave solution of (JHHJ). 



6 Numerical Example 

In the following, we specifically consider the third traveling solitary wave (see (fT|) 
solution of equation (jHHJ). 

Numerical values of the parameters are: Reh = 1.9, Ci = 3. 
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6.1 Analytical soliton 

The third traveling solitary wave solitary wave (see (JJJ) solution of equation (|55j). 
is given by: 

u(5, t) = Ux tanh [ C x (x - v t) ] + V (46) 

The variations of (}4l)j) a function of the non-dimensional space variable x and the 
non-dimensional time variable t for Ret = 1.9, Ci = 3, a = a 3 ~ 0.034, is displayed 
in Figure E 




Figure 1: The traveling solitary wave for Reh = 1.9, Ci = 3, a = 03 ~ 0.034 
Analytical calculation yield: 



1 

Re h 



2 C* 2 (8 C* 3 - 9 Ci )sech 2 [ d (x - v t)l r „ 
^ ; (8 C 2 - 13) 2 Re\ L 1 



L> 8 - 18i?e h + 2(8C| - 9Ci) tanh [Ci ( 

(47) 



As t increases, sup^ | ui + cuu^ — u^uix^ t) | decreases and tends towards 0; 
hence, the soliton (J4*6|) tends towards a solitary wave solution of (JIJ, with an error 

/ e -2 (2-u t) - 



i5| 



6.2 Finite difference calculation 

Advect n times the soliton (j4T)j) through the Finite second-order centered scheme 
in space, Euler-time scheme ETT1 in the space domain x £ [xq, +oo[, for: 

i. i?eh = 1.9, Ci = 3, cr = cr 3 ~ 0.034, which corresponds to the existence of a 
solitary wave solution of (jlj); 
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ii. Ren = 1.9, a = 0.06 = a ^ cr 3 , and Ren = 1.9, a = 0.07 = a" ^ cr 3 , which 
do not correspond to the existence of a solitary wave solution of (jlj). 

The mesh points number is equal to 150, with 100 points in the front wave. Xq is 

equal to —20 h. 

Denote: 

r = cr 3 , r = a , r = a (48) 

The numerical solitary wave, the analytical solitary wave (|46jh and the numeri- 
cal solutions at t = n t r = E [^^t 2 -] t = E [ ; ^rr] t , where £/ denotes the entire part, 
for n t = 0, n t = 5, n t = 10, n 4 = 50, respectively, are displayed in Fig. |21 




Figure 2: Numerical soliton, analytical soliton, and numerical solutions. Top: 
t = 0. Bottom: t = 10 r, t = 50 r. 
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